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Fingerprinting  Ammunition  on  a  Hybrid  Computer 


INTRODUCTION 

Computers  have  already  been  used  to  fingerprint  ammunition.  Geene^ 
states  that  his  model  provides  an  "aid  in  defining  more  meaningful  acceptance 
criteria  for  small  arms  propellants."  Geene's  study  was  aimed  at  illuminat¬ 
ing  the  reasons  for  ammunition  causing  weapon  performance  failures  after 
passing  acceptance  tests.  Geene's  "Energy  Equation"  model  does  provide  some 
insights  into  propellant  burn  and  projectile  dynamics,  however  his  model 
does  not  match  the  data  very  well. 

Instead  of  beginning  our  work  by  using  Geene's  model  as  a  springboard, 
our  approach  has  been  to  develop  a  set  of  non-linear  differential  equations 
which  describe  the  dynamics  of  the  physical  system.  The  simulation  of 
these  differential  equations  was  performed  on  an  EASE  2133  analog  computer 
and  a  PDP-7  digital  computer  connected  through  24  channels  of  analog/digital 
and  digital/analog  converters. 

SUMMARY 

This  report  summarizes  the  techniques  used  and  the  results  obtained  by 
fingerprinting  ammunition  on  a  hybrid  computer. 

A  FORTRAN  program  was  used  to  calculate  the  average  chamber  pressure 
versus  time  as  well  as  the  largest  and  smallest  pressure  used  for  the  averaqe 
of  ten  lots  of  ammunition  (20  rounds  per  lot).  The  data  points  were  taken 
from  200  photographs  furnished  by  Frankford  Arsenal.  Figures  la-le  show  the 
largest,  average  and  smallest  pressures  versus  time  for  each  of  the  furnished 
lots.  The  variation  between  the  largest  and  smallest  pressures  is  greater 
for  ball  propellant  (Lots  1,  2,  5-8,  10)  than  for  the  IMR  propellant 
(Lots  3,  4,  9). 
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A  hybrid  computer  model  was  developed  using  differential  equations  to 
describe  the  amount  of  propellant  burned,  the  pressure  buildup  and  decay, 
the  volume  behind  the  projectile,  and  the  acceleration,  velocity  and  posi¬ 
tion  of  the  projectile  within  the  muzzle.  The  coefficients  of  the  terms 
of  the  differential  equations  were  then  adjusted  so  that  the  model  curve 
matched  the  average  data  curves  for  the  ten  lots.  Table  1  shows  the  value 
of  the  coefficients  or  parameters  of  the  model  for  each  lot. 

Frankford  Arsenal  furnished  seven  unknown  pressure-time  curves  to  test 
the  fingerprinting  capability  of  our  approach.  Unfortunately ,  the  unknown 
data  did  not  come  from  the  same  lots  that  the  averages  were  computed  from. 
Nor  were  the  data  taken  at  the  same  ambient  temperatures.  Also,  it  would 
have  been  desirable  to  test  our  approach  more  thoroughly  through  many  more 
unknown  pressure-time  curves,  but  the  contract  time  expired.  The  results 
of  the  unknown  mathcings  are  shown  in  Figures  2a-2c. 

CONCLUSIONS 

The  fact  that  the  model  was  matched  to  each  lot,  as  verified  by  the 
results  shown  in  Table  1,  demonstrates  the  feasibility  of  our  approach. 

Although  the  unknowns  were  not  from  the  same  lots  as  furnished  origin¬ 
ally  by  Frankford  Arsenal,  some  interesting  observations  can  be  made. 

Figures  2a  and  2b  show  that  ambient  temperature  does  effect  the 
pressure-time  curve.  The  155°F  ambient  temperature  curve  of  Figure  2a  made 
a  very  good  match  with  Lot  9.  The  0°F  curve  fitted  Lot  9  best  but  not  nearly 
as  well  as  the  155  F  curve.  The  70°F  curve  fitted  Lot  3  better  than  Lot  9. 
Lots  3,  4,  and  9  of  the  original  data  were  IMR  tracer  lots  so  unknown 
TW-18179  is  an  IMR  propellant  lot. 
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TABLE  1  -  Equation  Coefficients 
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Figure  2b 

Unknown  Lot  WCC-6101 
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Figure  2c 
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The  curves  of  Figure  2b  show  that  at  -40°F  ambient  temperature,  the 
lot  matches  with  Lot  7,  but  at  155°F,  it  matches  with  Lot  9.  Lot  9,  as 
already  mentioned,  is  a  IMR  tracer  lot,  but  Lot  7  is  a  ball  propellant 
lot.  The  temperature  effect  of  unknown  WCC-6101  is  so  pronounced  that 
the  propellant  type  is  disguised. 

Figure  2c  shows  that  unknown  Lot  LC-12245  matches  with  Lot  5  which 
is  a  ball  propellant  lot  while  unknown  Lot  LC-12109  matches  with  Lot  9, 
an  IMR  tracer  lot.  Both  of  these  unknown  were  fired  at  70°F  ambient 
temperature.  All  of  the  original  lots  were  fired  at  70°F  ambient  temper¬ 
ature. 

Now  that  the  feasibility  of  this  approach  of- fingerprinting  the  chamber 
pressure  has  been  demonstrated  by  confirmation  that  the  unknown  lots  were 
correctly  identified,  it  would  seem  natural  to  extend  the  model  to  include 
the  effect  of  temperature,  and  chamber  pressure  on  the  port  pressure  and  then 
on  to  the  bolt  pressure.  A  total  model  would  then  show  the  direct  effects 
of  different  ammunition  on  bolt  action  and  the  important  parameters  could 
be  pinpointed. 

DISCUSSION 

Frankford  Arsenal  furnished  200  photographs  of  chamber  pressure  versus 
time.  There  were  20  photographs  for  each  of  ten  lots  of  ammunition.  The 
pressure-time  curve  was  divided  each  into  25  equal  segments  corresponding 
to  time  steps  of  0.05  milliseconds.  The  corresponding  pressure  for  each 
point  in  time  was  taken  from  the  photographs.  This  yielded  a  set  of  5,000 
data  points  which  were  punched  on  IBM  cards.  A  typical  pressure-time  curve 
as  furnished  by  Frankford  Arsenal  is  shown  in  Figure  3.  The  chopped  trace 
at  the  bottom  is  the  time  marker,  each  cycle  being  equal  to  0.1  milliseconds. 
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The  heavy  line  at  the  top  of  the  trace  is  the  60,000  psi  calibration  line 
for  the  chamber  pressure  curve  while  the  lower  heavy  line  is  the  20,000 
psi  calibration  line  for  the  port  pressure  trace. 

The  chamber  pressure  curve  is  characterized  by  an  initial  burning  of 
the  propellant  followed  by  a  rapid  increase  of  pressure  due  to  pressure, 
temperature,  burn-rate  relationships.  Part  way  up  the  pressure  trace, 
the  pressure  is  large  enough  to  overcome  the  equivalant  crimp  pressure  of 
the  cartridge  on  the  projectile  and  the  projectile  starts  to  move.  This 
causes  the  arc  of  the  curve  to  change  from  positive  to  negative  since  the 
volume  is  now  starting  to  increase.  The  pressure  is  still  rising  however 
and  will  until  the  volume  term  becomes  predominant.  The  initial  acceleration 
of  the  projectile  increases  with  pressure.  So  although  the  projectile  is 
still  accelerating  as  it  leaves  the  muzzle,  the  acceleration  is  less  as 
the  pressure  decreases . 

A  FORTRAN  program  was  written  which  calculated  the  average  pressure 
for  each  time  interval  as  well  as  the  largest  and  smallest  pressure  that 
went  into  the  average  calculation.  The  program  punched  IBM  cards  containing 
the  biggest,  smallest  and  average  pressures.  The  FORTRAN  program  was  run  on 
a  UN I VAC  1108  digital  computer. 

Figure  A1  of  the  Appendix  is  the  flow  chart  of  the  FORTRAN  program 
(Figure  A2).  The  IBM  cards  containing  the  780  pressure  data  points  are  read 
into  the  computer  and  immediately  printed  out  for  checking  purposes.  The 
average  for  each  point  in  time  is  computed  by  summing  the  pressure  points 
and  dividing  by  20.  As  the  average  is  being  computed  the  maximum  and  minimum 
pressures  that  went  into  the  computation  are  saved.  The  program  then  prints 
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a  list  of  the  largest,  smallest  and  average  pressure  for  each  point  in 
time.  Lastly,  a  set  of  IBM  cards  containing  the  largest,  smallest  and 
average  pressure  for  each  point  is  punched.  The  program  then  returns  to 
the  beginning  and  reads-in  the  next  set  of  data.  The  programs  continues 
through  this  loop  until  all  of  the  lots  have  been  inputed. 

The  FORTRAN  program  was  written  so  that  the  PDP-7  (digital  portion 
of  the  hybrid  computer)  would  not  have  to  store  5,200  data  points  in  core 
memory.  The  core  memory  contains  8,192  cells.  Instead,  the  data  has  been 
reduced  to  three  sets  of  26  points  per  lot  or  780  core  cells.  This  leaves 
a  larger  portion  of  memory  available  for  data  manipulation  and  model  simu¬ 
lation. 

For  a  complete  description  of  BNW's  Hybrid  Computer  I  and  our  software 
(SIMPL-1)  see  Reference  2.  That  paper  describes  our  approach  to  hybrid 
simulation. 

The  punched  cards  from  the  FORTRAN  program  were  read  into  the  hybrid 
computer  and  stored  in  data  tables.  By  selecting  switches  on  the  digital 
computer  console  any  one  of  the  ten  lots  may  be  displayed  along  with  the 
maximum  and  minimum  values. 

By  simply  using  the  averages  as  generated  by  the  FORTRAN  program,  the 
following  groupings  can  be  made  (see  Figures  la-le). 


Group 

I 

LOT# 

5,8 

Group 

II 

LOT# 

6,10 

Group 

III 

LOT# 

3,4,9 

Group 

IV 

LOT# 

1 

Group 

V 

LOT# 

2 

Group 

VI 

LOT# 

7 

The  average  of  Lots  5  and  8,  which  mate  up  Group  I,  are  almost  identical. 
The  two  curves  displayed  simultaneously  on  the  oscilloscope  almost  gives  the 
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appearance  of  being  one  curve.  There  exists  only  very  minute  differences 
between  the  averages.  Both  of  those  lots  are  ball  type  propellant. 

The  averages  of  Lots  6  and  10  which  make  up  Group  II  differ  in  the 
rate  of  increase  in  pressure  and  the  peak  pressure  attained.  They  agree  in 
that  the  peak  pressure  occurs  at  the  same  time  and  the  pressure  decay  for 
both  is  practically  identical.  Lot  6  is  a  tracer  with  ball  propellant  and 
Lot  10  is  a  ball  propellant. 

The  averages  of  Lots  3,  4,  and  9  which  compose  Group  III,  are  grouped 
together  because  they  are  very  similar  in  rate  of  rise  in  pressure,  peak 
pressure,  time  of  peak  pressure  and  pressure  decay.  Although  these  lots  are 
very  similar  each  one  has  a  different  rate  of  rise  in  pressure  and  a  different 
peak  pressure.  They  are  however  so  similar  in  appearance  that  all  three  can 
be  within  the  Large  -  Small  limits  of  any  one  of  them.  Lots  3,  4  and  9  are 
tracers  with  IMR  propellant. 

The  averages  of  Lots  1 ,  2  and  7  which  compose  Group  IV,  Group  V,  and 
Group  VI  respectively,  are  not  directly  similar  to  any  other  lot  previously 
grouped  or  with  any  other  lot  mentioned  here.  Each  one  is  a  case  of  its  own. 
Lot  1  could  be  grouped  with  Lots  5  and  8  since  the  curves  are  very  nearly 
identical  except  that  Lot  1  is  delayed  by  .05  milliseconds.  Lot  7  could  be 
grouped  with  Lot  5  and  8  since  they  all  have  very  similar  pressure  decay 
curves  and  initial  pressure  rise  curves.  Lot  7  however  has  a  higher  rate 
of  increase  in  pressure  and  reaches  a  higher  peak  pressure  because  the  time 
of  peak  pressure  for  Lot  5  and  8  and  Lot  7  are  identical.  Lot  2  is  truly 
an  independent.  This  lot  is  not  close  to  any  other  group.  Its  rate  of 
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increase  in  pressure  is  much  slower,  the  peak  pressure  occurs  much  later, 
and  the  peak  pressure  is  higher  than  any  other  lot.  Lot  1 ,  2  and  7  are 
ball  propellant  lots. 

Model  Equations 

The  current  model  has  been  constructed  from  simple  physical  relation¬ 
ships  and  some  generalized  assumptions. 

B  is  defined  as  the  volume  of  propellant  bed  burned  and  has  units  of 
cubic  inches.  The  rate  of  change  of  B  is  described  by 

6  =  kjT  +  k2P 

where  B  is  the  derivative  of  B  with  respect  to  time  or  dB/dt.  6  increases 
as  T  and  P  increase.  If  the  assumption  is  made  that 

T  =  k] "b 


(1) 


(2) 


then  equation  (1)  can  be  written 


where 


B  =  k-jB  +  k2P 
kl  =  klkl 


(3) 


In  reality  the  proportionality  constant  k-j  relates  the  amount  of 
propellant  burnt  to  the  temperature  which  is  then  related  to  the  burn- 
rate.  In  the  model  k^  relates  the  amount  of  propellant  burnt  to  the  burn- 
rate  and  has  the  dimension  sec"^ .  The  proportionality  constant  k2  relates 
the  pressure  to  the  burn-rate  and  has  dimensions  of  in5/lb-sec.  This  is 
similar  to  the  a  used  by  Geene^  in  his  burn-rate  equation. 

P  is  the  chamber  pressure  in  psig  and  is  defined  by 

p  =  k3  V  +  k4B  +  P0 


or 


p  =  k3  7  f  '  k3  ' f  k4fi 


(4) 
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In  this  model  the  pressure  is  derived  from  two  sources.  One  source 
is  the  ideal  gas  law  that  relates  the  pressure  of  the  existing  gases  to 
their  temperature  and  volume.  The  other  is  due  to  the  production  of  the 
gases  in  the  combustion  of  the  propellant.  If  the  substitutions 

II 

T  =  k1  B 
f  =  k/ft 

I  II 

k3  =  k3  k1 
\l  =  Ax 


The  term  Pq  is  an  initial  pressure  due  the  primer.  The  proportional¬ 
ity  constant  k3  is  usually  seen  in  the  form  of  the  ideal  gas  law  nR  where  n 
is  the  number  of  moles  of  gas  and  R  is  the  ideal  gas  law  constant.  In  the 
model  k3  has  the  dimension  lb/in  .  The  proportionality  constant  k^  relates 
the  amount  of  propellant  burnt  to  pressure.  In  fact,  it  might  be  considered 
similar  to  an  efficiency  term  because  it  says  that  for  so  much  burnt  pro¬ 
pellant  the  pressure  should  increase  by  so  much.  In  the  model  it  has  the 
dimension  of  lb/in  , 

The  equations  describing  the  volume  behind  the  projectile  are: 
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where  VT  is  the  total  volume  behind  the  projectile,  VQ  is  the  initial 
voids  volume  within  the  cartridge,  B  is  the  volume  of  propellant  burned 
and  is  the  volume  of  the  muzzle.  The  units  of  volume  are  cubic  inches. 
The  volume  of  the  cartridge  was  given  as  0.113  in3  and  the  specific  volume 
of  the  propellants  as  65-70%.  Using  70%  as  the  specific  volume  of  the 
propellant  Vq  =  (.11 3) (.3)  =  0.034  in3  and  therefore  the  maximum  the  B 
can  be  is  0.079  in3. 

The  volume  of  the  muzzle  can  be  described  by  =  0.038  x  where 
0.038  is  the  cross-sectional  area  of  the  muzzle  in  square  inches,  x  is 
the  length  of  the  muzzle  in  inches.  The  maximum  value  for  x  is  18.1  inches; 
x  is  zero  until  the  projectile  starts  to  move,  x  which  is  a  function  of 
time  can  be  obtained  by  the  following  equations 


where  x  is  the  acceleration  of  the  projectile  in  inches/sec2  and  k  is 

d 

a  parameter  converting  the  units  of  pressure  to  units  of  acceleration. 

can  also  be  used  to  obtain  the  measured  muzzle  velocity  (approximately 
38,400  in/sec)  since 


The  units  of  k  are  in/lb-sec2.  P  is  the  chamber  pressure  in  psig  and  P 

sw 

is  the  pressure  at  the  time  the  projectile  starts  moving.  P$w  is  a 
parameter  which  can  be  varied  in  the  model. 
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To  obtain  x,  x  is  integrated 


t 


where  kv  is  a  psuedo-friction  term  and  is  unitless,  k  can  be  used  to 

X  X 

adjust  x  to  be  18.1  inches  at  the  appropriate  time,  k  and  k  are  para- 

A  a 

meters  which  can  be  varied  in  the  model. 

The  error  is  determined  by 

t  9 

c  f  (AVRAGE  -  P)  c  -j 

E  -  /  — p-~ap—  dT 

0 

where  AVRAGE  (digital  computer  mnemonic)  is  the  average  data  curve  for 
any  of  the  ten  lots.  The  AP  term  in  the  denominator  is  small  but  necessary 
It  is  used  to  keep  the  division  from  blowing  up  at  t  =  0  when 
(AVRAGE  -  P)/P  =  0/0.  The  model  is  matched  to  AVERAGE  when  E  is  minimized. 

The  analog  computer  diagram  as  well  as  the  digital  program,  digital 
flow  diagram  and  digital  memory  allocation  is  shown  in  the  Appendix 
(Figures  A3-A6).  The  circles  are  symbolic  of  analog  potentiometers.  Most 
of  these  are  available  on  the  analog  computer  console  and  they  are  used  for 
rapid  adjustment  of  the  parameters.  The  triangles  with  a  box  on  one  of 
the  edges  are  analog  integrators.  The  hexagons  are  trunk  lines  which 
connect  the  analog  signals  (MX's)  to  the  digital  computer  and  the  digital 
signals  (DA's)  to  the  analog  computer. 

An  advantage  of  analog  computers  is  that  they  integrate  differential 
equations  exactly  and  in  parallel.  One  of  the  advantages  of  the  digital 
computer  is  that  it  accurately  performs  arithmetic  functions. 


The  fingerprinting  model  is  set  up  to  use  these  advantages.  The 
digital  computer  Is  programmed  to  solve  all  of  the  differential  equations, 
display  the  results  and  change  parameter  values  rapidly. 

Model  Constants 

As  the  model  now  stands  these  constants  are  only  an  indication  of 
some  real  physical  constant.  In  general,  they  give  an  indication  of  the 
magnitude  of  the  effect  that  a  parameter  has  upon  the  actual  physical 
system.  The  actual  number  and  its  units  have  no  real  physical  significance, 
for  most  of  the  constants,  at  this  time. 

The  term  Pq  is  used  by  the  model  to  get  started.  It  is  related  to  the 
pressure  output  of  the  primer.  It  however  should  not  be  considered  as  the 
actual  pressure  output  of  the  primer,  but  rather  the  pressure  due  to  the 
primer  as  seen  by  the  sensing  pressure  gauge.  This  gauge  is  some  distance 
away  from  the  primer  and  much  of  the  initial  force  of  the  primer  would  be 
absorbed  by  the  propellant  between  the  primer  and  the  gauge.  The  pressure 
output  of  the  primer  as  reported  by  Squire  and  Devine^  was  measured  in  a 

t 

closed  bomb  at  approximately  30,000  psig.  The  model  begins  with  a  pressure 
of  approximately  350  psig. 

The  constant  termed  switching  pressure  Pgw  or  the  pressure  at  which 
the  motion  of  the  projectile  becomes  a  significant  factor  in  the  pressure- 
time  relationship,  appears  to  have  much  more  physical  significance.  Squire 

(3) 

and  Devine'  '  report  that  a  valid  estimate  of  this  pressure  is  18,000  psig. 
The  model  predicts  pressures  ranging  from  13,000  to  25,000  psig  for  this 
pressure.  Lot  2  showed  a  switching  pressure  of  31,000  psig.  This  is  not 
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out  of  order  because  Lot  2  has  a  much  later  and  higher  peak  pressure  than 
any  other  lot.  Also,  if  one  looks  carefully  at  the  photo  for  Lot  2  there  is 
an  inflection  point  in  the  data  at  30,000  psig  where  the  slope  of  the 
pressure  curve  decreases,  which  is  what  would  happen  as  the  volume  began  to 
increase  as  the  projectile  moves.  If  one  looks  at  the  switching  pressure 
used  by  the  model  and  then  the  data,  there  is  this  noticeable  inflection 
point  in  almost  all  of  the  lots. 


The  proportionality  constants  kx  and  ka  are  primarily  used  to  calibrate 
the  model.  They  are  to  be  adjusted  such  that  the  distance  parameter  x  would 
be  about  18  inches.  The  length  of  the  barrel  is  the  only  available  data 
for  the  calibration  of  these  constants.  Data  concerning  the  acceleration 
or  velocity  of  the  projectile  at  any  time  while  it  is  in  the  muzzle  was  not 
available.  For  this  first  model  too  many  simplifications  had  been  made  and 
the  lack  of  data  meant  that  this  calibration  could  not  be  accomplished.  If 
the  source  of  term  x  is  written  in  this  fashion,  one  can  see  the  possibility 
of  combining  the  constants  k  and  k  into  one  constant. 

X  a 


This  is  pointed  out  also  by  the  fact  that  if  the  multiplication  of  k  ,  k  , 

X  fit 

and  P  are  performed  and  integrated  by  the  model,  x  has  the  following  form 

x  =  k  P  t2  .14  <  k  <  .16 

For  all  practical  purposes  this  proportionality  constant  k  is  constant. 
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The  proportionality  constant  k.j  was  held  constant  while  the  data  for 
Table  1  was  taken.  This  was  done  because  the  model  was  relatively  insensi¬ 
tive  to  k-|.  The  response  of  the  model  to  a  large  change  in  k-j  could  be 
completely  masked  by  a  small  change  in  k2.  What  this  says  is  that  the  burn- 
rate  is  relatively  insensitive  to  temperature  once  combustion  has  started. 

The  way  that  the  model  currently  exists,  the  proportionality  constants 
k2j  k3»  k4  interact  a  great  deal.  This  interaction  makes  it  impossible  to 
directly  separate  the  physical  significance  of  any  one  given  constant. 
Currently  the  constants  k|,  k2,  kg,  k4  are  not  unique  for  any  given  lot 
number.  If  the  coefficients  k^  and  k2  could  be  set  according  to  some  cali¬ 
bration  which  says  that  so  much  propel laht  will  be  burnt  at  some  time,  this 
would  then  remove  some  of  the  freedom  of  the  constants  and  tend  to  produce 
unique  values  for  k  -j ,  k2,  kg,  k^.  These  constants  would  then  have  more  direct 
physical  significance.  Data  wasn't  available  to  accomplish  this.  The  model 
was  sensitive  to  the  proportionality  constant  kg.  This  is  the  constant  that 
turns  the  pressure  around  and  relates  the  volume  to  the  pressure.  Of  the 
constants  k^ ,  k2,  kg,  k^,  the  constant  kg  is  the  most  independent  constant 
of  the  group.  If  kg  is  used  to  separate  the  10  lots  into  groups  the  following 
grouping  could  be  made 

Group  A  Lots  1,2,6,7,10 

Group  B  Lots  3,  4,  9 

Group  C  Lots  5,8 

Group  B  is  the  same  as  Group  III,  and  Group  C  is  the  same  as  Group  I  as  pre¬ 
viously  defined  by  the  photographs  of  the  average  data. 

The  constants  k2  and  k^  both  control  the  upward  swing  of  the  pressure. 

Due  to  magnitudes  in  the  model  the  positive  slope  portion  of  the  pressure 
curve  is  almost  completely  controlled  by  the  product  of  k2  k4  P.  This  means 
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that  if  kg  is  increased  then  k^  should  reflect  some  equivalent  decrease. 

This  is  evident  by  comparing  the  two  very  similar  lots,  Lots  6  and  10. 

For  Lot  6,  kg  is  3.231x10“^  and  k^  is  1 .386x10^  while  for  Lot  10,  kg  is 
-6  5 

4.393x10  and  k^  is  1.178x10  .  This  same  type  of  thing  is  evident  in 
Lots  3,  4,  9. 

Model  Modifications 

The  current  simplified  model  only  partially  relates  directly  to  the 
real  physical  system.  With  a  few  modifications  one  could  still  have  a  rela¬ 
tively  simplified  model  yet  which  would  relate  more  directly  to  the  physical 
system. 

One  such  modification  would  be  to  include  temperature  T  in  the  model 
explicitly  rather  than  implicitly.  This  equation  would  be  a  boundary  value 
equation  of  the  form 

t  =  -kT 

T(0)  =  0  (11) 

TH  =  Ta 

Here  T  is  the  adiabatic  flame  temperature  of  the  propellant. 

a 

An  additional  modification  can  be  made  in  equation  (3)  for  the  pressure. 

In  equation  (3)  the  ideal  gas  law  written  as  k 3  y-  •  The  constant  kg'  is 
usually  seen  as  nR  where  n  is  the  number  of  moles  of  gas  and  R  is  the  ideal 
gas  constant.  Since  the  number  of  moles  of  gas  changes  as  the  propellant 
is  burnt,  terms  in  the  model  need  to  reflect  this  increase. 

When  the  temperature  T  and  moles  of  gas  n  are  included  in  the  model, 
the  model  would  then  become 

e  =  k,  t  +  k2  p 

P  =  y-  (nf  +  Til  -  y-V)  +  k^  6  (12) 

t  ■  -  k6  T 
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Some  more  major  type  of  modifications  would  include  dividing  the 

propellant  bed  into  sections  of  different  characteristics.  These  sections  . 

(3) 

are  described  by  Squire  and  Devine v  .  Such  a  model  would  include  the 
differences  in  propellant  burning  throughout  the  propellant  bed,  the  effect 
of  compaction  on  the  burn-rate,  the  geometry  of  the  burning  face  of  the 
propellant,  the  effect  of  pressure  differentials  on  the  burning  rate,  and 
the  effect  of  volume  change  on  temperature,  and  burn-rate.  Another  parameter 
that  needs  to  be  considered  is  the  muzzle  temperature.  Such  a  model  would 
indeed  be  much  more  complex  and  require  more  time  to  implement. 

Since  the  model  up  to  this  point  has  started  with  conditions  after  the 
primer  has  done  its  job,  perhaps  some  modeling  should  be  done  from  the  time 
the  firing  pin  strikes  the  primer  to  the  point  where  the  existing  model  begins 
This  type  of  model  would  predict  the  condition  of  the  propellant  bed,  tempera¬ 
ture,  pressure  and  other  parameters  due  to  the  primer.  By  combining  this 
model  with  the  previous  one,  the  model  should  be  such  that  it  could  relate  to 
the  actual  physical  system  from  the  time  the  firing  pin  strikes  the  primer  to 
the  time  the  projectile  leaves  the  muzzle. 
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Figure  A2 
FORTRAN  Program 
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Figure  A3 

Analog  Computer  Diagram 
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Figure  A4 

Digital  Computer  Program 
(Hybrid) 
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PROGRAM  NAME  :  FINGER 
PAGE  1  OF  LISTING 


FINGER 

/  LEA  1-12-70 

/  MAIN  PROGRAM  LOOP 

BGN,  JMS  INIT  /  INITIALIZATION  PROGRAM 

JMS  CLEAR  /  CLEARS  ALL  D  TO  A  LINES 
LOOP,  LAW  DKMTBL  /  DKM$  INSTRUCTION  TABLE 

JMS  DKM.  /  DKM$ 

JMS  CHECK  /  EXAMINES  AC  SWITCHES 

JMS  GETMOD  /  GET  OPERATION  MODE 

SNA  /  POTSET? 

JMP  LOOP-1  /YES 

Al,  JMS  DISPLA 

JMS  MODEL 
LAW  INTI 
JMS  RINT1 . 

0 

ISZ  4777 
JMP  LOOP 
JMP  LOOP 

RELOC=5000  /  RELOCATION  CONSTANT 


INTI,  SX  0 


PROGRAM  NAME  :  FINGER 
PAGE  2  OF  LISTING 


/  TA3LE  OF  CHARACTERS  AND  LOCATIONS  FOR  DKM$  CONTROL 

DKMTBL,  220  /  CONTROL  P 

RESTRT  /  RE-I  NITIALIZE  MAIN  PROGRAM 

224  /  CONTROL  T 

LOADDT  /  LOAD  DDT 

HLT 

/  PROGRAM  TO  LOAD  DDT 

LOADDT.  DPI  /  DISABLE  PRIORITY  INTERRUPT 

MI  CL  /  CLEAR  MODIFIED  INSTRUCTIONS 

LAC  (12000  /  UNIT  1 

MMSE  /  SELECT  TRANSPORT 

LAC  (-’10020 
IACISZA 

JMP  .-1  /  35  MS  DELAY 

LAC  (76  /  BLOCK  NUMBER 

DAC  17701 

JMS  1  7620  /  READ  FILE 

MMLC  10  /  STOP  TAPE 

LAC  BTMSYM  /  BOTTOM  OF  SYMBOL  TABLE 

DAC  15216 

LAC  POISV  /  POINTER  LOCATION 

DAC  15204  /  RESET  POINTER 

JMP  16000  /  START  DDT 
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PROGRAM  NAME  :  FINGER 
PAGE  3  OF  LISTING 


/  PROGRAM  FOR  CLEARING  ALL  D  TO  A  LINES 
CLEAR,  0 


CLA 

DA0 

/  CLEAR  AC 
DAI 

DA2 

DA3 

DA4 

DA5 

DA6 

DA  7 

DA  8 

DA9 

DAI  3 

DAI  1 

DA12 

DAI  3 

DAI  4 

DAI  5 

DAI  6 

DAI  7 

DAI  8 

DAI  9 

DA2  0 

DA21 

DA22 

DA23 

JMP  I 

CLEAR 

/INITIALIZATION  ROUTINE 

I  NIT,  2 

LAC  (-’13000 
I  AC!  SZA 

JMP  .-1  /  DELAY  BEFORE  CAF 

CAF  /  CLEAR  ALL  FLAGS 

LAC  15216  /  ADDRESS  OF  BOTTOM  OF  DDT’S  SYMBOL  TABLE 

DAC  BTMSY# M  /  SAVE 

LAC  1 5204  /  DDT  POINTER  LOCATION 

DAC  POISV#  /  SAVE  DDT  POINTER  LOCATION 

LAC  (20000  /  UNIT  2 

MMSE  /  SELECT  TRANSPORT 

LAC  (-’10000 
I  AC! SZA 

JMP  .-1  /  35  MS  DELAY 

LAC  (2  /  PROGRAM  DDTFIL  (ALL  SUBROUTINES) 

DAC  1  7701 

LAC  (RELOC  /  RELOCATION  CONSTANT 
DAC  17702 

JMS  1  7620  /  READ  FILE 

MMLC  10  /  STOP  TAPE 

LAC  BTMSYM 

DAC  SYMBL.  /  SYMBLS 
LAC  POISV 

DAC  POINT.  /  POINTS 
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PROGRAM 
PAGE  A 

RESTRT, 


NAME  *  FINGER 
OF  LISTING 


CAF 


CAC 

MIST 

JMS  GETMOD 
XOR  (4 
SZA 

LAC  (3 
JMS  SETMOD 


/  MODIFIED  INSTRUCTIONS 
/  GET  OPERATION  MODE 

/  IS  IT  COMPUTE? 

/YES  SET  TO  IC 
/  SET  OPERATION  MODE 


LAM 

DAC  #SWMNTR 
DAC  #MXCNTR 
LAW  DATA 
DAC  DATADR# 
LAW  KONTBL 
DAC  #CONTBL 
DZM  TT# 


/  SET  COUNTER  TO  CHECK  AC  SWITCHES 


/  ADDRESS  OF  TABLES 

/  THIS  IS  THE  TABLE  WHERE  THE  MODEL  CONSTANTS 
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MXRN 

/  RANDOM  MULTIPLEXING 

JMS  ICL. 

/  ICLS 

TS, 

1 

/  TIME  SCALE 

FC , 

5000* 

/  CLOCK  FREQUENCY 

NIC, 

50* 

/  #  OF  ITERATIONS  IN  COMPUTE 

NIIC, 

50’ 

/  #  OF  ITERATIONS  IN  IC 

CLSCMP 

/  CLOCK  SERVICE  COMPUTE 

CLSIC 

/  CLOCK  SERVICE  IC 

0 

/  HOLD  (NO  ROUTINE) 

TIM, 

0 

/  TIME  SINCE  ENTERING  COMPUTE 

LAC  (4000  /  ADOV 

ASC  /  ENABLE  A  TO  D  OVERLOAD 

EPI  /  ENABLE  PRIORITY  INTERRUPT 


LAW  TBUF  /  TELETYPE  BUFFER  TABLE  ADDRESS 
LMQ 

LAW  TBUFE  /  END  OF  TELETYPE  BUFFER 

JMS  IATYP .  /  INITIALIZE  PRIORITY  TYPING  (IATYP$> 

JMP  I  INIT 

TBUF -RELOC  3542 
TBUF E=RELOC  3717 

/  PROGRAM  FOR  SETTING  OPERATION  MODE 
SETMOD,  0 

DAC  MODE.  /  SAVE  MODE  (SMODE) 

LAC  MODFLG 
XOR  (ANAL 

SZA  /  IS  MODE  ANALOG  OR  DIGITAL? 

JMP  I  SETMOD  /  DIGITAL 

LAC  MODE.  /  GET  MODE  (SMODE) 

JMS  SMODE.  /  SET  ANALOG  MODE  (SMODES) 

JMP  I  SETMOD 

/  PROGRAM  FOR  DETERMINING  OPERATION  MODE 

GETMOD,  0 

LAC  MODFLG 
XOR  (ANAL 

SNA  /  IS  MODE  ANALOG  OR  DIGITAL? 

JMS  AMODE.  /  ANALOG  (AMODES) 

LAC  MODE.  /  DIGITAL  (GET  SMODE) 

JMP  I  GETMOD 

MODFLG,  ANAL  /  DIGT  OR  ANAL 

DIGT=12345 

ANAL=01234 
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/  VARIABLES  USED  IN  VARIABLE  STORING  ROUTINES 

VARX,  1  /  NUMBER  OF  VARIABLES  TO  BE  STORED 

VAR1,  TIME 

VAR2,  TIM 

VAR3,  TIME 

VAR  4,  TIME 

VAR  5,  TIME 

VAR6,  TIME 

/  CLOCK  SERVICE  ROUTINES 

CLSCMP, 

21/  HLT  /  FOR  ILLEGAL  CAL 

44/  NOP  /  A  TO  D  OVERLOAD 

CLSCMP/  LAS 

AND  (4 000  /  LOOK  AT  SWITCH  6 

YflR  #1  STSW 

SNA  /  TIME  TO  STORE  VARIABLES? 

JMP  CLSIC  3'  /NO  USE  IC  ROUTINE 

LAC  ENT1# 

SMA  /  TIME  =  ZERO? 

JMP  .+4  /  NO 

LAC  TIM  /  YES 

DAC  TREF  /  STORE  REFERENCE  TIME 

DZM  ENT1  /  SET  FLAG 

LAC  TIM 

SUB  TREF#  /  SUBTRACT  REFERENCE 

DAC  TIME  /  GIVES  VALUE  OF  TIME 

STL  /  SET  LINK  FOR  STORES 

JMS  STORE.  /  STORES 

TAU1  /  LOCATION  OF  TABLE  OF  TIME  INTERVALS 

STRT3L  /  STORAGE  TABLE  ADDRESS 

NSTOR,  300’  /  #  OF  LOCATIONS  IN  TABLE 

LAC  NIC  /  #  OF  ITERATIONS  PER  SECOND  (COMPUTE  MODE) 
JMP  RCL . 

CLSIC,  LAS 

AND  (4000 
DAC  LSTSW 
LAM 

DAC  ENT l 
DZM  TIME# 

CLL 

JMP  NSTOR-3 
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/  TIME  INCREMENT  TABLE  FOR  STORAGE  PROGRAM  STORES 


TAU1 , 

'10 

TIMEl  , 

10*00 

TAU2, 

1  *00 

TIME2 , 

50*00 

TAU3, 

5*00 

TIME3, 

1000*00 

TAU4, 

*01 

TIME4 , 

1100*00 

TAU5 , 

0 

TIME5, 

0 

TAU6, 

0 

TIMES, 

0 

TAU7, 

0 

TIME7, 

0 

.  10V 


/  INTERVAL  FOR  STORAGE  UNTIL  TIMEl 


/  STORAGE  TABLE  FOR  STORES 


STRTBL  ,  0 

.  300'/ 
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/  AC  SWITCH  CHECKING  PROGRAM 

CHECK,  0 

LAS 

AND  (60000 

SNA  /  TIME  TO  SET  MODE?  (SWITCHES  3  AND  4) 

JMP  VRCHK  /  NO 

CAS  (40000  /  YES 

JMP  .+  6  /  HOLD 

JMP  .+3  /  IC 

LAC  (4  /  COMPUTE 

JMP  .+  4 

LAC  (3 

SKP 

LAC  (1 

JMS  SETMOD  /  SET  OPERATION  MODE 
VRCHK,  LAS 

CLL 

AND  (1000  /  LOOK  AT  SWITCH  8 

SZA  /  TIME  TO  LIST  TABLE? 

STL  /  YES 

LAW  STRTBL  /  NO 

JMS  VRTYP .  /  VRTYPS 

NTYP,  300’ 

CLL 

LAS 

AND  (2000  /  LOOK  AT  SWITCH  7 

SZA  /  TIME  TO  PRINT  VARIABLES  CONTINUOUSLY? 

STL  /  YES 

JMS  VMNTR.  /  VMNTRS 
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ISZ  CNT123  /  COUNTER  FOR  SETTING  VARIABLE  LOCATIONS 
JMP  I  CHECK 
LAC  VARX 

DAC  VARX.  /  VARXS 
LAC  VAR1 

OAC  VAR1  .  /  VAR1  $ 

LAC  VAR2 

DAC  VAR2  .  /  VAR2S 

LAC  VAR3 

DAC  VAR3 .  /  VAR3S 

LAC  VARA 

DAC  VARA.  /  VARA$ 

LAC  VAR 5 

DAC  VAR5 .  /  VAR5S 

LAC  VAR6 

DAC  VAR6.  /  VAR6$ 

LAM- 700’ 

DAC  CNT123 
JMP  I  CHECK 

CNT123,  LAM-20* 
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ICL.=REL0C  2601 
DKM.=RELOC  562  1 
VRTYP.=RELOC  1 72  1 
STORE. =RELOC  21  1 
VMNTR.=R ELOC  432  1 
TRIP1 .  =R ELOC  2146  1 
FG2.=REL0C  2457 
SMODE.=R  ELOC  2545 
AMODE.=RELOC  2523 
HMPY1  .=R  ELOC  2266 
HMPY2.=REL0C  2316 
HDIV.=RELOC  2400 
X .  1  .  =R ELOC  2444 
TFRM.=RELOC  3360 
TMSG  ,=R ELOC  3235 
HMPY3.=REL0C  2347 
RINTl.=RELOC  1  777  1 
VAR1  .  =RELOC  553  l 
VARX.  =  VAR1  .-1 
VAR2.  =  VAR1  .+  1 
VAR3.  =  VAR2.+  1 
VAR4.  =  VAR3.+  1 
VAR5.  =  VAR4.+  1 
VAR6.  =  VAR5.  + 1 
TTAB.=RELOC  3367 
TCR.=RELOC  3347 
TVOLT.=RELOC  3156 
HEADR.=RELOC  502  l 
TTYP.=RELOC  3535 
TFLG ,=R ELOC  3520 
TOCT.=RELOC  3053 
TTAG.sRELOC  1677  l 
ATYP  »=R  ELOC  3424 
TDEC.=RELOC  3101 
POINTER  ELOC  1676  l 
SYMBL.=R  ELOC  1675  1 
TSP.=RELOC  3340 
CL2.=REL0C  3051 
XI  0.  =RELOC  2424 
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END.^RELOC  2127  1 
FC.=RELOC  3046 
TS.=RELOC  3047 
MODE*  =R ELOC  3052 
CL1 » =RELOC  3050 
RCL.=RELOC  3021 
IATYP.=RELOC  3376 

NCP^NOP 
SNA  =SNA 

hlt=hlt 

CMA =CMA 
SMA=SMA 
SPAzSPA 
SZA=SZA 
SKP=SKP 
CLL=CLL 
3TL=STL 
CLArCLA 
LAS=LAS 

lmq=lmq 

LA  CQ  =  LA  CQ 
LRSS-LRSS 
LLSS=LLSS 


PROGRAM  NAME  S  FINGER 
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/  THIS  IS  THE  WORKING  PORTION  OF  FINGER 
/  THIS  IS  THE  DISPLAY  AND  THE  MODEL 


/  DISPLA  IS  A  FUNCTION  GENERATOR  PROGRAM  THAT  DISPALYS  THE  FOLLOWING 


DATA  BY  CONTROL  OF  THE  AC  SWITCHES 

LARGEST,  SMALLEST,  AVERAGE,  FOR  ANY  GIVEN  LOT  OF  DATA 
THE  AC  SWITCHES  SELECT  THE  LOT  BEING  VIEWED 


/AC  SWITCHES 


LOTS 


10  DECIMAL  OR 


OCTAL 


DISPLA, 


LAC  TT 
CAS  Cl  .0 
JMP  DOFUNC 
NOP 

ISZ  FLAG# 
ISZ  SWMNTR 
JMP  DOFUNC 
LAM  -20’ 
DAC  SWMNTR 
LAC  LOTNUM 


/  >  SO  IN  THE  MIDDLE  OF  A  CALCULATION 
/  = 

/  <  SO  SET  FLAG  AND  LOOK 

/  THIS  DECIDES  TO  SEE  IF  WANT  DIFFERENT  LOT 
/  DON’T  CHECHK  NOW  JUST  DO 
/  RESET  COUNTER 

/  NEED  TO  SEE  IF  ANYONE  IS  IN  THERE 
/  IF  0  THEN  FIRST  TIME  THROUGH 

/  CHECK  TO  SEE  IF  YOU  WANT  TO  CHANGE  LOTS 
/  BIT  13  IS  THE  ONE  1  YES,  0  NO 


DOFUNC 


ADD  (1 
CAS  LOTNUM 
SKP 

JMP  DOFUNC 
DAC  LOTNUM 
CLL 
MUL 
116 
LACQ 

ADD  (-116 
ADD  DATADR 
DAC  #ADRESD 
LAC  (-114 
DAC  COUNT# 
LAW  TABLAD 
DAC  #ADREST 
LAC  I  ADRESD 
DAC  I  ADREST 
ISZ  ADRESD 
ISZ  ADREST 
ISZ  COUNT 
JMP  ,-5 


/IF  YOU  FORGET  YOU  GET  LOT  1 
/  DO  YOU  ALREDY  HAVE  THIS  ONE  LOADED 

/  YEP  NO  NEED  TO  RELOAD  A  NEW  LOT 
/  THIS  IS  TO  CHECK  NEXT  TIME 


/  THIS  IS  THE  LOCATION  OF  THE  START  OF  DATA 
/  THIS  IS  THE  POINTER  USED  IN  LOADING 

/  THIS  KEEPS  TRACK  OF  WHEN  TO  STOP  LOADING 

/  THIS  IS  TH  POINTER  OF  WHERE  TO  LOAD 


/  DONE 
/  NOPE 
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DOFUNC,  MX0  /  GET!  THE  TIME 

JMS  ADFLG 
DAC  TT 
•  LAM  -4+1 
JMS  FGMODS 
LAW  TABLE 
LAC  TT 

DA0  /  SEND  OVER  THE  LARGEST 

DAI  /  SEND  OVER  THE  SMALLEST 

DAC  #A VRAGE  /STORE  THE  AVERAGE  FOR  THE  FUTURE 

DAC  #UNKWN  /  THIS  IS  THE  UNKNOWN 

DA6 

LAC  AVRAGE 

DA2  /  SEND  OVER  THE  AVERAGE 

JMP  I  DISPLA 


TABLE,  0 

2500. 
LRS  6 
25’ 

TABLAD,  0 

TABLAD  78*/ 
UNKADD,  0 

UNKADD  28V 


ADFLG,  0  /  THIS  IS  MX  LINE  CHECK 

ADSF 
JMP  .-1 
ADRB 

JMP  I  ADFLG 

OVERLD,  0  /  THIS  STORES  ADDRES  WHERE  DIGITAL  OVERLOAD 

722102 
LAC  OVERLD 
DAC  OVERAD 
JMP  I  OVERLD 
OVERAD,  0 
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MODEL,  0  /  THIS  IS  THE  MODEL  PORTION  OF  THE  PROGRAM 

LAC  FLAG  /  FLAG=1  IF  START  OF  NEW  CALCULATION 
SNA  /  FLAG=0  IN  THE  PROCESS  OF  CALLCULATING 

JMP  DOMODL 
DZM  FLAG 

ISZ  MXCNTR  /  ONLY  SEE  IF  CONSTANTS  ARE  HOW  ONCE  IN  A 
JMP  DOMODL  /  WHILE  NOT  EVERY  TIME 
LAM-5’  /  THIS  TIME  CHECK 

DAC  MXCNTR 

LAS  /  BIT  10=1  THEN  GET  THE  ONES  IN  MEMORY 

AND  (200 

SNA 

JMP  WRTCON 

LAC  LOTNUM  /  THIS  IS  THE  LOT  NUMBER  OF  DATA  WANTED 

SUB  (1 

CLL 

MUL 

10 

LACQ 

ADD  CONTBL 
DAC  ADRESD 
LAW  MODCON+ 1 
DAC  ADREST 
LAC  (-7 
DAC  COUNT 
LAC  I  ADRESD 
DAC  I  ADREST 
ISZ  ADREST 
ISZ  ADRESD 
ISZ  COUNT 
JMP  .-5 
JMP  DOMODL 
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WRTCON,  LAS  /  THIS  IS  TO  SEE  IF  YOU  WANT  TO  WRITE 

AND  (100  /  BIT  11=1  THEN  STORE  THESE  AWAY  FOR  THIS  LOT 

SNA 

JMP  GETCON 
LAC  LOTNUM 

HLT  /  TO  MAKE  SURE  YOU  ARE  WRITTING  AND  KNOW  IT 

LAS  /IF  WANT  TO  WRITE  MUST  HAVE  LOTNUM  ON  AC  SW 

AND  <17 
CMA 

AND  LOTNUM 
SZA 

JMP  DNTWRT  /  NOT  THE  SAME  SO  DO  NOT  WRITE 
LAC  LOTNUM 
SUB  (1 
CLL 
MUL 
10 

LACQ 

ADD  CONTBL 
DAC  ADREST 
LAW  MODCON+1 
DAC  ADRESD 
LAC  (-7 
DAC  COUNT 
LAC  I  ADRESD 
DAC  I  ADREST 
ISZ  ADREST 
ISZ  ADRESD 
ISZ  COUNT 
JMP  .-5 
JMP  DOMODL 

DNTWRT,  LAC  LOTNUM  /  THIS  IS  FOR  NOT  WRITING 
CMA 

HLT  /  TO  LET  YOU  KNOW  THAT  NO  WRITE 

JMP  DOMODL 
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GETCON, 

KK1  , 

KK2, 

KK3, 

KK4, 

KP0, 


LAS 

AND  (40  /  BIT  12=1  THEN  MX 

S  NA 

JMP  DOMODL  /  NO  READ ,  WRITE, MX  SO  USE  WHAT  IS  THERE 

MX1  /  RELATES  BDOT  ,B 

J MS  ADFLG 
DAC  K1 

MX2  /  RELATES  BDOT,P 

JMS  ADFLG  ' 

DAC  K2 

MX3  /RELATES  T,P,B 

JMS  ADFLG 
DAC  K3 


MX4  /  RELATES  BDOT ,PDOT 

JMS  ADFLG 
DAC  K4 

MX5  /  THE  INITIAL  PRESSURE 

JMS  ADFLG 
DAC  P0 


KKX, 


MX6 

JMS  ADFLG 
DAC.  KX 


/  PSEUDO  FRICTION  CONSTANT 


KPSW, 


MX  7 

JMS  ADFLG 
DAC  PSW 


/  THE  SWITCHING  PRESSURE 


KKA, 


MX8 

JMS  ADFLG 

DAC  KA  /  THE  ACCERERATION  CONSTANT 


JMP  DOMODL  /  THIS  IS  TO  GET  OVER  THE  CONSTANT  TABLE 
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MODCON,  0 

K1 ,  0 

K2,  0 

K3,  0 

K4,  0 

P0,  0 

KX,  0 

PSW,  0 

KA,  0 


/  THIS  IS  DUMMY  SO  TABLE  WILL  HAVE  TAG 
/  1 0f 4  SF 
/  1 0 1 7  SF 
/  .5  SF 
/  .0005  SF 
/  .1  SF 
/  10f5  SF 
/  PSW  .001  SF 
/  1/SEC  10t4  SF 


/  LIST  OF  OTHER  SCALE  FACTORS 
/  DIGITAL  ANALOG 

/  b  INT3  2000  SF  1000  SF 

/  BOOT  I N t3/S£C  1000  SF  1000  SF 

/  p  LB/I Nt2  .001  SF 

/  PDOT  LB/I NT2 /SEC  .001  SF 

/  v  INt3  2000  SF 

/A  I  Nt2  10f4  SF 

/  X2D0T  IN/SECT2  .04  SF 

/  XDOT  IN/SEC  .04  SF 

/  X  IN  5  SF 


/  NOW  WITH  CONSTANTS  LETS  DO  SOME  WORK 


DOMODL,  MX 9  /  GET  B  THE  AMOUNT  OF  MATERIAL  BURNT 

JMS  ADFLG  /  1000  B 
ALSS+1 

DAC  B#  /  2000B 
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OF  LISTING 

CBDOT, 

LAC 

P 

/ 

.001  P 

JMS 

X10. 

/ 

.01P 

SZL 

JMS 

OVERLD 

JMS 

HMPY2  . 

LAC 

K2 

/ 

10T7  K2 

SZL 

JMS 

OVERLD 

DAC 

SPOT 1 # 

/ 

1 000 (K2*P) 

LAC 

B 

/ 

2000  B 

JMS 

HMPY2 • 

LAC 

K 1 

/ 

1 0 14  K 1 

SZL 

JMS 

OVERLD 

JMS 

HMPY1  . 

/ 

2*1 0t5  K1*B 

2.5 

/ 

1000  K 1*B 

ADD 

SPOT1 

/ 

1000(K1*B+K2*P) 

DAC 

BOOT# 

/ 

1000  *  BDOT 

JMS 

X.l . 

/ 

100  BDOT 

CMA 

/ 

-100  BDOT 

DAC 

# DA BOOT 

CX2  DOT, 

LAC 

P 

/ 

.001  P 

JMS 

HMPY2. 

LAC 

KA 

/ 

1 0 f 4  KA 

SZL 

/ 

.1  K  A*P 

JMS 

OVERLD 

JMS 

HMPY1  . 

40. 

DAC 

X2  DOT 

/ 

.04  X2D0T 

56 


PROGRAM  NAME  t  FINGER 
PAGE  19  OF  LISTING 


CPSW,  LAC  TT  /  THIS  IS  SUCH  ONLY  SWITCH  XDOT  ON  NOT  OFF 

CAS  <0.1  /  WHEN  IN  COMPUTE  s  ,  '  '  v*r 

NOP  /  >  0.1  '  .:-:r 

JMP  .+  6  /  =  0.1  N. ,  -v 

DZM  #S WITCH  /<  0.1 

DZM  XDOT  .  c.  ' . v.  *. 

DZM  X  ^  \  . 

DZM  P  . 

JMP  .+  10’  -  s? 


CV, 


LAC 

SWITCH 

"  r  v  ?  • 

CAS 

(0 

SKP 

/ 

>0  '  'K:V  :  > 

SKP 

/ 

=  0  '  , 

JMP 

.+8’ 

/ 

<0 

LAC 

P 

r  /\  \  i  )fi  y  wjk  s-  .*  • 

CAS 

PSW 

V  \  •  ;  V 

JMP 

.+4 

/ 

>  psw 

NOP 

/ 

=  PSW 

DZM 

X2  DOT 

/ 

<  PSW 

SKP 

ISZ 

SWITCH 

*  U  X  '«■  A  H  "  T  0  C  S  }  u  it  £  ia » 

?oae  fiSFij 

LAC 

XDOT 

/ 

.04  XDOT  y  ; 

JMS 

HMPY2  • 

LAC 

KX 

/ 

12t5  sf,93, iSi> 

SZL 

JMS 

OVERLD 

/ 

-40*XDOT*KX 

mi 

mih.i 

' ,  \ 

;  'J 3/ >* A  *  C--  T  00  0  )  A  X  ;  f,  \  , 

12.  t 

?p»  . 

SZL 

'>*  '  I  j'  ,  • 

JMS 

OVERLD 

•  i  %y  5  ■(  ’0; 

DAC 

KXXDOT 

/ 

5  XDOT*KX 

LAC  X  /  5X 

ALSS+2  /  20  X 

JMS  HMPY1  . 

380.  /  MUZZLE  AREA  <  INt2  )  -10t4  SF 

SZL  '  -  • 

JMS  OVERLD 

ADD  (68.  /  INITIAL  CHAMBER  VOLUME  2*10t3  SF 

ADD  B  /  BURNT  MATERIAL  (  INT3  2*l0f3  SF 

DAC  V#  /  2000  V 
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CP  DO  T ,  LAC  XDOT  /  .04  XDOT 

JMS  HMPY1. 

380.  /  10f4  A 

SZL 

JMS  OVERLD  /  4*A*XD0T 
JMS  HMPY3. 

LAC  B  /  2000  B 

LAC  V  /  2000  V 

SZL 

JMS  OVERLD 

DAC  SP0T2#  /  4*B*A*XD0T/V 
LAC  BDOT  /  1000  BDOT 

ALSS+2  /  4000  BDOT 

JMS  X.l .  /  400  BDOT 

JMS  X.l.  /  40  BDOT 

JMS  X.l.  /  4*BD0T 

SUB  SP0T2 

JMS  HMPY3  .  /  4*<BD0T-B*A*XD0T/V> 

LAC  K3  /  K3/2 

LAC  V  /  2000V 

SZL 

JMS  OVERLD 

DAC  SP0T2  /  .001  K3 (BD0T-8*A*XD0T/V)/V 

LAC  BDOT  /  1000  BDOT 

JMS  HMPY3  • 

LAC  K  4  /  .0005  K4 

LAC  (500. 

SZL  /  .001  K4*BD0T 

JMS  OVERLD 
ADD  SP0T2 

DAC  SP0T2  /  .001K3 (BDOT-B*A*XDOT/V)/V+.001K4*BDOT 

LAC  P0  /  .1  P0 

JMS  X.l.  /  .01  P0 

JMS  X.l .  /  .001  P0 

ADD  SP0T2 

DAC  PDOT  /  .001  PDOT 

LAC  P 

NOP 

DA4 


LAC  DABDOT  /  -100  BDOT 
DA3 
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C  ER  R  , 

ERR, 


SX, 

X2DCT, 

KXXDOT , 

PDOT, 


NAME  :  FINGER 
OF  LISTING 


LAC  AVRAGE  /  THIS  IS  TO  GET  RID  OF  THE  0  IN  THE  DENOM 

JMS  ERR 

DA5 

JMP  I  MODEL 
0 

DAC  N#  UMRAP 
ADD  DELP 
DAC  #DENOMP 

LAC  NUMRAP  /  START  GETTING  THE  ERROR. 

SUB  UNKWN 
JMS  HDIV. 

LAC  DENOMP 
SZL 

JMS  OVERLD 
JMS  HSQ$ 

SZL 

JMS  OVERLD 
JMS  X.l  . 

JMS  X.l. 

CMA 

JMP  I  ERR 
JMP  I  MODEL 

0.  /  IC  FOR  X2D0T 

0. 

LRSS+2 
ADD  XDOT# 

0 

DAC  XDOT 

0.  /  IC  FOR  XDOT 

0. 

LRSS+0 
ADD  X# 

0 

DAC  X 


0.  /  IC  FOR  P 

0. 

LRSS+0 
ADD  P# 

0 

JMP  .+1 
DAC  P 
JMP  END. 
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DATA, 


0.00 

0.50 

2.00 

3.50 

5.00 

5.50 

7.30 

10.00 

15.00 

23.00 

32.00 

40.50 

47.00 

51  .00 

52.00 

51  .00 

49.00 

45.00 

39.50 

35.50 

31.00 

27.00 

23.50 

20.00 

18.00 

16.00 

0.00 

0.00 

0.00 

0.00 

,  0.00 

0.00 

1  .03 

2.33 

5.50 

11.00 

19.50 

23.00 

33.52 

36.00 

40.00 

44.00 

42.00 

36.50 

31  .33 

26.00 

22.00 

19.00 

17.00 

14.50 

12.50 

11.50 

0.33 

2.05 

1 .05 

1 .72 

2 . 70 

4.03 

5.52 

7.55 

11.57 

17.95 

25.52 

32.67 

39 . 45 

44.30 

48.20 

48.47 

45.52 

40.92 

35.35 

32.55 

26.12 

22.70 

19.55 

17.23 

15.15 

13.47 

3.33 

2.00 

2.00 

3.53 

5.50 

6.50 

8.3'3 

10.50 

1  6.00 

21 .50 

27.00 

32.50 

39.50 

43.02 

56.00 

62.00 

62.50 

61 .00 

56.00 

43.32 

41  .00 

35.50 

31 .00 

27.00 

24.00 

21  .23 

2.33 

0.30 

0,00 

1 .00 

2.50 

4,50 

5.50 

6.53 

7,50 

10.00 

13.50 

18.50 

23 . 50 

23.50 

33.50 

38.00 

43.00 

46.00 

42 .02 

35.50 

30.50 

26.00 

22.50 

19.50 

17.53 

1  5.33 

0.02 

3.23 

3.72 

2.00 

3.72 

5.67 

6.60 

8.32 

10.55 

15.00 

20.85 

27.03 

33 . 1  5 

42.00 

46.37 

51  .45 

53.92 

51  .72 

47.42 

43.37 

35.15 

30.30 

25  ,95 

22.42 

19.65 

1  7.27 

3.00 

0.50 

2.00 

3.00 

8.30 

12.50 

22.52 

32.50 

42.00 

49.50 

53.00 

54.00 

53.33 

48.23 

42.00 

36.00 

31  .50 

27.00 

23.00 

23.00 

13.00 

15.50 

14.00 

12.50 

11.50 

12.50 

0.32 

0.00 

3.00 

0.53 

1  .50 

3.00 

10.03 

19.50 

31  .50 

38.50 

45.50 

47.50 

44.30 

40.00 

34.50 

30.00 

25  .5  0 

22.00 

18.00 

1  6.00 

14.00 

12.00 

10.50 

9.50 

8.50 

3.00 

0.00 

0.30 

3.95 

2.03 

4.32 

7.77 

15.92 

27.35 

37.27 

44.15 

49.10 

50.52 

48.08 

43.65 

37.92 

32.55 

28.00 

24.05 

23.45 

17.83 

15.52 

13.55 

1  1.97 

10.65 

9.52 

8.60 

3.32 

0.50 

1  .50 

4.00 

9.50 

12.50 

22.50 

33.00 

41  .50 

48.00 

52.00 

52.50 

49 .50 

46.50 

43  .50 

.  39.50 

35  .00 

30.00 

26.50 

23.00 

20.00 

17.00 

15.00 

13.50 

12.32 

10.50 
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0.20 

0.00 

3.20 

1.00 

2.50 

6.00 

10.32 

13.00 

25.00 

31.50 

39.00 

43.00 

44.32 

37.22 

32.00 

28.00 

24.00 

20.50 

17.52 

15.00 

13.50 

12.00 

10.50 

9.50 

8.50 

8.00 

3.32 

0.03 

0.60 

2.10 

5.60 

9.65 

15.52 

25.35 

33.50 

40.12 

45  .85 

48.47 

47.77 

44.37 

39.67 

34.57 

30.05 

25.70 

22.15 

19.07 

16.55 

14.47 

12.70 

11.32 

12.37 

9.12 

3.30 

0.00 

1  .50 

3.00 

6.00 

9.50 

12.52 

21.50 

32.00 

37.50 

43.00 

48.00 

52.00 

52.23 

52.00  ■ 

51.00 

48.00 

42.50 

36.53 

31.50 

2  6.50 

22.50 

19.50 

17.50 

15.50 

14.00 

3.20 

0.00 

0.00 

0.50 

1 .00 

3.00 

5.50 

7.30 

10.00 

16.00 

24.00 

30.50 

37.50 

44.50 

45  .00  • 

40.00 

34.50 

30.00 

26.00 

22.00 

19.50 

l  7.00 

14.50 

13.5  0 

12.00 

10.50 

3.00 

2.00 

0.62 

1  .67 

3.13 

5.78 

8.55 

12.97 

19.95 

27.22 

34.50 

40.00 

44.95 

43.05 

48.05 

44.75 

39.95 

34.85 

29.77 

26.03 

22.30 

19.35 

16.97 

15.07 

13.40 

12.02 

0.00 

2.00 

1  .00 

2.50 

7.50 

12.00 

22.00 

33.50 

41  .50 

49.30 

54.00 

58.00 

59.50 

56.00 

54.00 

47.50 

42.00 

35.50 

30.30 

26.00 

22.50 

19.50 

17.00 

15.00 

13.53 

12.30 

0.00 

0.00 

0.00 

0.50 

1.50 

2.50 

5.50 

9.50 

16.50  . 

26.00 

34.50 

42.00 

46.00 

46.00 

38.50 

32.00 

28.00 

23.50 

20.00 

17.50 

15.50 

13.50 

12.00 

10.50 

9.50 

8.50 

0.00 

0.00 

0.30 

1.35 

3.05 

6.07 

10.82 

19.30 

29.02 

38.02 

45.35 

50.72 

52.80 

52.92 

46.02 

39.62 

34.08 

28.80 

24.72 

21  .23 

18.50 

16.27 

14.27 

,  12.67 

11.32 

10.22 

0.00 

1  .00 

1 .50 

3.50 

5.00 

7.50 

11.50 

21  .00 

31  .50 

41  .00 

48.50 

53.00 

55.50 

54.00 

52.00 

48.00 

42.00 

36.50 

31  .03 

26.50 

23.00 

20.00 

18.00 

16.00 

14.00 

13.00 

3.33 

0.00 

0.00 

0.50 

1 .50 

2.00 

3.00 

7.00 

13.50 

24.00 

32.50 

40.00 

44.50 

47.50 

45.00 

39.00 

34.00 

28.50 

24.50 

20.50 

13.00 

16.00 

14.00 

12.00 

11  .03 

10.00 
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0.00 

0.30 

0.92 

1.72 

3.15 

4.80 

7.37 

11 .72 

20.85 

30.45 

38.62 

45.67 

49.55 

51  .22 

48.80 

43.80 

37.95 

32.70 

27.83 

23.85 

20.53 

18.00 

15,80 

14.02 

12.55 

11.25 

1.03 

2.00 

3.00 

5.00 

6.50 

9.00 

13.50 

22.50 

30.00 

38.50 

44.50 

49.00 

51  .00 

51  .50 

50.00 

47.00 

42.50 

37.50 

32.00 

27.53 

24  .00 

21.00 

20.00 

16.00 

14.53 

13.00 

-1  .33 

-1  .00 

0.00 

0.50 

1.50 

2.50 

4.50 

8.50 

14  .50 

22.50 

29.00 

35.50 

40.00 

44.00 

44  .50 

38.50 

33.50 

28.50 

24.53 

21.50 

19.00 

16.50 

15  .00 

13.50 

12.00 

11.00 

0.37 

0.40 

1  .35 

2., 2  5 

3.80 

5.50 

7.72 

11.82 

18.87 

26.52 

33.95 

40.45 

45.15 

47.70 

47.55 

44.02 

39.45 

34.00 

29.20 

25.20 

21 .80 

18.70 

16.87 

14.80 

13.17 

.  11.80 

0.00 

0.00 

1 .50 

4.50 

8.00 

14.00 

24.53 

35.50 

43.00 

48.50 

51.50 

51 .00 

49.30 

45.00 

41  .00 

36.50 

32.00 

28.00 

24.50 

21  .50 

19.00 

16.50 

14.50 

13.00 

11.50 

10.50 

0.03 

0.00 

0.50 

2.00 

5.00 

8.00 

12.00 

23.00 

32.00 

39.00 

44.00 

47.00 

45.50 

40  .03 

35.00 

31  .00 

26.50 

23.00 

20.03 

17.00 

15.00 

13.00 

11.50 

10.50 

9.50 

8.00 

0.00 

0.00 

0.90 

2.95 

6.72 

10.35 

20.23 

30.70 

39.12 

45.55 

49.12 

49.52 

47.30 

42.95 

37.97 

33.35 

28.88 

24.92 

21 .65 

18.77 

16.38 

14.35 

12.72 

11.40 

10.17 

9.12 

3.00 

2.00 

4.00 

7.00 

11.50 

19.00 

27.50 

34.50 

42.50 

49.00 

54.00 

55.00 

53.50 

53.00 

51.00 

47.50 

42.50 

36.50 

31.00 

27.00 

23.00 

20.00 

17.50 

15.50 

14.00 

12.50 

0.00 

0.00 

0.00 

0.50 

1 .00 

2.00 

5.00 

10.00 

17.00 

23.50 

31.00 

38.50 

44.50 

44.00 

38.00 

32.00 

28.00 

24.00 

21.03 

18.00 

15.50 

13.50 

12.00 

10.50 

9.50 

8.00 

0.00 

0.65 

1  .67 

3.07 

5.32 

8.42 

13.87 

20.92 

28.90 

37.00 

43.27 

48.47 

50.87 

49.25 

44.90 

38.90 

33.80 

28.82 

24.72 

21  .42 

18.60 

16.35 

14.60 

12.87 

1 1 . 62 

10.47 
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/  KONTBL  IS  THE  TABLE  OF  CONSTANTS  FOT  THE  MODEL 

KONTBL,  0 

KONTBL  80V 

LOTNUM,  0 

DELP,  0.5 

PAUSE 


63 


FI-NGER 

ADFLG 

1315=*= 

ADRESD 

3473=*= 

ADREST 

3474* 

AMODE. 

7523 

ANAL 

1234 

ATYP. 

10424 

AVRAGE 

3475* 

A1 

10* 

B 

3476* 

3D0T 

3477=*= 

BGN 

0* 

BTMSYM 

3500* 

CBDOT 

1  505* 

GERR 

1663* 

CHECK 

766* 

CLEAR 

47* 

CL  S  CMP 

224* 

CLSIC 

251* 

CL1  . 

10050 

CL2 . 

10051 

CNT123 

1045* 

CONTBL 

3501* 

COUNT 

3502* 

CPDOT 

1612* 

CPSW 

1542* 

cv 

1567* 

CX2  DOT 

1532* 

DABDOT 

3503* 

DATA 

1735* 

DATADR 

3504* 

DELP 

3472* 

DENOMP 

3505* 

DIGT 

12345 

DISPLA 

1046* 

DKMTBL 

22* 

DKM. 

5563 

DNTWRT 

1427* 

DOFUNC 

1120* 

DOMODL 

1501* 

END. 

7130 

ENT1 

3506* 

ERR 

1  667* 

FC 

155* 

FC. 

10046 

FGMODS 

FG2 . 

7457 

FLAG 

3507* 

GETCON 

1433* 

GETMOD 

205* 

HD  I V  • 

7400 

HEADR. 

5503 

HMPY1  . 

7266 

HMPY2 . 

7316 

HMPY3  . 

734  7 

HSQS 

IATYP  • 

10376 

ICL. 

7601 

INIT 

102* 

INTI 

20* 

KA 

1500* 

KKA 

1464* 

KKX 

1456* 

KK1 

1437* 

KK2 

1442* 

KK3 

1445* 

KK4 

1450* 

KONTBL 

3351* 

KPSW 

1461* 

KP0 

1453* 

KX 

14  76* 

KXXDOT 

1720* 

K 1 

1471* 

K2 

1472* 

K3 

1473* 

K4 

1474* 

LOADDT 

27* 

LOOP 

2* 

LOTNUM 

3471* 

LSTSW 

3510* 

MODCON 

1470* 

MODEL 

1330* 

MODE. 

10052 

MODFLG 

214* 

MXCNTR 

3511* 

NIC 

156* 

NIIC 

15  7* 

NSTOR 

246* 

NTYP 

1013* 

NUMRAP 

3512* 

OVERAD 

1327* 

OVERLD 

1322* 

P 

3513* 

PDOT 

1726* 

POINT. 

6677 

POISV 

3514* 

PSW 

1477* 

P0 

1475* 

RCL. 

10021 

RELOC 

5000 

RESTRT 

132* 

RINTl  . 

7000 

SETMOD 

1  74* 

SMODE. 

7545 

SPOT1 

3515* 

SP0T2 

3516* 

STORE. 

5022 

STRTBL 

311* 

SWITCH 

3517* 

SWMNTR 

3520* 

SX 

1711* 

SYMBL. 

6676 

TABLAD 

1143* 

TABLE 

1137* 

TAU1 

261* 

TAU2 

263* 

TAU3 

265* 

TAU4 

267* 

TAU5 

271* 

TAU6 

273* 

TAU7 

275* 

TBUF 

10542 

TBUFE 

10717 

TCR. 

10347 

TDEC. 

10101 

TFLG  • 

10520 

TFRM. 

10360 

TIM 

163* 

TIME 

3521* 

TIME1 

262* 

TIME2 

264* 

TIME3 

266* 

TIME4 

270* 

TIME5 

272* 

TIME6 

274* 

TIME7 

276* 

TMSG, 

10235 

TOCT. 

10053 

TREF 

3522* 

TRIP1 . 

7147 

TS 

154* 

TSP. 

10340 

TS. 

10047 

TT 

3523* 

TTAB. 

10367 
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FI'NGER 

TTAG  . 

67i30 

TTYP. 

1  0535 

TVOLT • 

10156 

UNKADD 

1261  * 

UNKWN 

3524* 

V 

3525* 

VARX 

21 5^ 

VARX. 

5553 

VAR1 

216* 

VAR1. 

5554 

VAR2 

21 7* 

VAR2  • 

5555 

VAR3 

220* 

VAR3  • 

5556 

VAR  4 

221* 

VAR4.  . 

5557 

VAR5 

222* 

VAR  5. 

5560 

VAR  6 

223* 

VARS. 

5561 

VMNTR . 

5433 

VRCHK 

1  004* 

VRTYP. 

5173 

WRTCON 

1370* 

X 

352  6* 

XDOT 

3527* 

X.l  . 

7444 

XI 0. 

7424 

X2D0T 

1712* 

FI-NGER 


ANAL 

1234 

TBUF 

10542 

CNT123 

1045* 

RELOC 

3000 

TBUFE 

10717 

DISPLA 

1046* 

STORE. 

5022 

DIGT 

12345 

DOFUNC 

1120* 

VRTYP . 

5173 

BGN 

0* 

TABLE 

113  7* 

VMNTR. 

5433 

LOOP 

2* 

TABLAD 

1143* 

HEADR. 

5503 

A1 

10* 

UNKADD 

1261* 

VARX  . 

5553 

INTI 

20* 

ADFLG 

1315* 

VAR1  . 

5554 

DKMTBL 

22* 

OVERLD 

1322* 

VAR2 . 

5555 

LOADDT 

27* 

OVERAD 

1327* 

VAR3  • 

5556 

CLEAR 

47* 

MODEL 

1330* 

VAR  4. 

5557 

INIT 

102* 

WRTCON 

1370* 

VAR  5. 

5560 

RESTRT 

132* 

DNTWRT 

1427* 

VAR6. 

5561 

TS 

154* 

GETCON 

1433* 

DKM . 

55  63 

FC 

155* 

KK1 

143  7* 

SYMBL  • 

6676 

NIC 

156* 

KK2 

1442* 

POINT. 

6677 

NIIC 

157* 

KK3 

1445* 

TTAG  . 

6700 

TIM 

1 63* 

KK4 

1450* 

RINT1  . 

7000 

SETMOD 

174* 

KP0 

1453* 

END. 

7130 

GETMOD 

205* 

KKX 

1456* 

TRIP1  . 

7147 

MODFLG 

214* 

KPSW 

1461* 

HMPY1  . 

72  66 

VARX 

215* 

KKA 

1464* 

HMPY2  . 

7316 

VAR1 

216* 

MODCON 

1470* 

HMPY3. 

7347 

VAR2 

217* 

K  1 

1471* 

HDI  V. 

7400 

VAR3 

220* 

K2 

1472* 

X10. 

7424 

VAR4 

221* 

K3 

1473* 

X  .1 . 

7444 

VAR5 

222* 

K4 

1474* 

FG2. 

7457 

VAR6 

223* 

P0 

1475* 

A  MODE. 

7523 

CLSCMP 

224* 

KX 

1476* 

SMODE. 

7545 

NSTOR 

246* 

PSW 

1477* 

I  CL. 

7601 

CLSIC 

251* 

KA 

1500* 

RCL. 

10021 

TAU1 

261* 

DOMODL 

1501* 

FC. 

10046 

TIME1 

262* 

CBDOT 

1505* 

TS. 

10047 

TAU2 

263* 

CX2DOT 

1532* 

CL1  . 

10050 

TIME2 

264* 

CPSW 

1542* 

CL2. 

10051 

TAU3 

265* 

CV 

1567* 

MODE. 

10052 

TIME3 

266* 

CPDOT 

1612* 

TOCT. 

10053 

TAU4 

267* 

CERR 

1663* 

TDEC . 

10101 

TIME4 

270* 

ERR 

1 667* 

TVOLT . 

101  56 

TAU5 

271* 

sx 

1711* 

TMSG  . 

10235 

TIMES 

272* 

X2D0T 

1712* 

TSP  • 

10340 

TAU6 

273* 

KXXDOT 

1720* 

TCR . 

10347 

TIME6 

274* 

PDOT 

1726* 

TFRM . 

10360 

TAU7 

275* 

DATA 

1735* 

TTAB . 

103  67 

TIME7 

276* 

KONTBL 

3351* 

IATYP. 

10376 

STRTBL 

311* 

LOTNUM 

3471* 

ATYP. 

10424 

CHECK 

766* 

DELP 

3472* 

TFLG  . 

10520 

VRCHK 

1004* 

ADRESD 

3473* 

TTYP. 

10535 

NTYP 

1013* 

ADREST 

3474* 
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FI-NG  ER 

avrage 

3475* 

B 

3476* 

BDOT 

3477* 

STMSYM 

3500* 

CONTBL 

3501* 

COUNT 

3502* 

DABDOT 

3503* 

DA  TADR 

3504* 

DENOMP 

3505* 

ENT  1 

350  6* 

FLAG 

3507* 

LSTSW 

351  0* 

MXCNTR 

3511* 

NUMRAP 

3512* 

P 

3513* 

POISV 

3514* 

SPOT1 

3515* 

SP0T2 

351  6* 

SWITCH 

351  7* 

SWPINTR 

3520* 

TIME 

3521* 

TREF 

3522* 

TT 

3523* 

UNKWN 

3524* 

V 

3525* 

X 

3526* 

XDOT 

HSQ$ 

FGMODS 

352  7* 
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Figure  A5 

Digital  Computer 
Flow  Diagram 
(Hybrid) 
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Calculate  the  address  of  the 
data  for  this  Lot  Number 


Load  the  data  for  this  Lot  Number  into  the 
dynamic  buffer  used  by  the  Arbitrary  Function! 
Generator 


Use  the  Arbitrary  Function  Generator  to 
Display  the  data  for  this  Lot  Number  and 
the  unknown 


£ 


Calculate  the  address  of 
the  model  constants  for 
this  Lot  Number 


Load  the  model  constants 
in  the  dynamic  buffer 
into  the  storage  table 


Calculate  the  address  of 
the  model  constants  for 
this  Lot  Number 


Load  these  constants  into 
the  dynamic  buffer  for 
constants  used  in  the  model! 


ever  been 
>  P 

v  sw^- 

X  =  0. 


No^ 

P  >  P 

sw 

Yes  . 

DA  the  value  ERRM 


a 
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Figure  A6 


Digital  Computer 
Memory  Map 
(Hybrid) 
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Absolute  Memory  Addresses 


3600g  h. 


3500.  -► 
8 


Digital  Computer  Data  Storage  Map 

Model  Constants  Lot  #  10 

Model  Constants  Lot  #  9 

Model  Constants  Lot  #  8 

Model  Constants  Lot  #  7 

Model  Constants  Lot  #  6 

Model  Constants  Lot  #  5 

Model  Constants  Lot  #  4 

Model  Constants  Lot  #  3 

Model  Constants  Lot  #  2 

Model  Constants  Lot  #  1 

.Average 

.Small  Lot  #  10 

Large 

"  ■"  1  ' '  1  >  . " -  . 

_Average 

.Small  Lot  #  9 

Large _  ’ _ 

^Average 

.Small  Lot  #  8 

Large _ _ _ 

.Average 

.Small  Lot#  7 

Large _ _ _ 

Average 

Small  Lot  #  6 

Large  _ _ _ _ 

Average 

Small  Lot#  5 

Large _ 

Average 

Small  Lot  #  4 

Large _ _ 


